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Phase Structure of Lattice QCD at Finite Density with Dynamical 
Fermions 

I.M. Barbour a , J.B. Kogut b and S.E. Morrison" * 

a Dept. of Physics and Astronomy, University of Glasgow G12 8QQ, U.K., (UKQCD Collaboration) 
b Dept. of Physics, University of Illinois, 1110 West Green Street, Urbana, IL 61801 

We compare the chemical potential associated with the onset of non-zero baryon number density on 6 4 and 
8 4 lattices at (3 = 5.1 and ma=0.01. We provide evidence for Z(3) tunnelling. We determine a critical chemical 
potential of {xa ~ 0.1 which is unexpectedly low. We also determine the dependence of the onset of the observed 
phase transition on the quark mass. The physically misleading result of the quenched theory is shown to persist 
despite the inclusion of the complex fermion determinant. 



1. Introduction 

Lattice QCD at finite baryon density holds the 
key to an understanding of the phase transition 
between quark-confined hadronic matter and the 
quark-gluon plasma. At zero temperature we ex- 
pect this transition to occur at /j, = corre- 
sponding to the lowest lying state with non-zero 
baryon number. 

The fact that the fermion determinant becomes 
complex when [i is non-zero makes numerical sim- 
ulations particularly difficult. Simulations using 
the quenched theory (|] give the physically mis- 
leading result that as the quark mass tends to 
zero, the lowest lying state with non-zero baryon 
number becomes massless. This corresponds to 
fj, c ~ Recent work by Stephanov |!| using 

the random matrix method has provided a solu- 
tion to this particular problem by demonstrating 
that at nonzero /x the quenched theory is not a 
simple n — * limit of a theory with n quarks. 
This explanation implies that the phase of the 
determinant is important for a simulation of full 
QCD at finite /i. However earlier work on the full 
theory at strong coupling || indicates that the 
transition still occurs at 

1.1. Formulation 

The GCPF can be expressed as an ensemble 
average of the fermionic determinant normalised 



with respect to the fermionic determinant at fi = 
0. This means that our probability weighting fac- 
tor is well-defined. 



_ $[dU][dW}detM(ii)e- s ^ u > u ^ 
~ j[dU][dW]detM{^ = 0) e SAU,u^] 

which leads to: 
/ detM (/i) 



Z = 



\ detM (/j, = 0) 



(1) 



(2) 



where S g is the pure gauge action and M is the 
Kogut-Susskind fermion matrix: 



iM = G + im + Ve/ 1 + Vh^. 



(3) 



G contains all the spacelike gauge links and the 
bare quark mass and V all the forward timelike 
links. On each given configuration of the gauge 
fields detM is expanded as a polynomial in the 
fugacity variable with powers ranging between 
[— 3n^,3n^] (n s corresponding to the spatial ex- 
tent of the lattice). The ensemble averages of 
the coefficients of this normalised expansion are 
the corresponding expansion coefficients for the 
GCPF. 

The gauge field configurations are generated 
using Hybrid Monte Carlo, with the expansion 
coefficients of detM obtained from the eigenval- 
ues of the propagator matrix P (Gibbs B) 
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Thus (for a = 1): 

det(iM) = e 3 ' tn ' nt det(P-e-' i ). (5) 

det(P-e-»)= w n e~ nli (6) 

The eigenvalues, A, of P have a Z(nt) symme- 
try and a A and 1/A* symmetry, which leads to 
an expansion of the form: 

3n s 3 3n 3 3 

Z = < b H> e nAmt = e- (£ "-™^ /T (7) 

n— — 3n s 3 n— — 3n s 3 

1.2. Z(3) Tunnelling 

In the pure gauge theory we have N c equivalent 
vacua related by Z(N C ) rotations. Tunnelling be- 
tween the different Z(3) vacua is much more prob- 
able in the confined sector than in the deconfined 
sector. Since the pure-gauge action as well as the 
integration measure are invariant under the Z(3) 
transformation, the GCPF can also be written as: 

, . = J[dU}[dW}detM(fi + z 3 /n t )e- s °V' u ^ 
W J[dU}[dm}detM(n = 0)e- s ^ u > u ^ () 

Averaging over the three Z(3) vacua would elim- 
inate the triality non-zero coefficients and, since 
Z(0) = 1, the sum of the triality zero coefficients 
should tend to 1. As Fig. 1 shows, at ft = 5.1 and 
m = 0.01, we obtained a clear signal for Z(3) tun- 
nelling. The sum of the triality zero coefficients 
shows a strong tendency to average to 1 whilst 
those of the triality one and two coefficients each 
tend to zero. This is consistent with simulating 
in the confined phase at fi = as expected at this 
coupling and quark mass. 

1.3. Baryon Number Density 

We consider a description of the system in 
terms of the canonical partition functions for 
fixed particle number. This technique allows us to 
include only those triality zero coefficients which, 
at the end of the simulation, are positive and have 
a reasonable error. The chemical potential as a 
function of the baryon number density is obtained 
from the local derivative of the energy, e n , of the 
state with n fermions with respect to n. 

"<"» = sj£ < 9 » 



where p is the fermion density, 

In Fig. 2 we show the number density, at 
(3 = 5.1, m = 0.01 on 6 4 and 8 4 lattices, as a 
function of chemical potential. Only the depen- 
dence at low densities is shown as it is only the 
coefficients < b n > for n small which are deter- 
mined with reasonable statistical error. There is 
little difference between the 6 4 and the 8 4 number 
densities. The onset of non-zero baryon number 
density at /i c — 0.1 is unexpectedly low. 

1.4. Mass Dependence of the Onset fi 

Therefore we have also considered the chemi- 
cal potential (the onset fi) required to make the 
level with 3 quarks equally probable with the zero 
quark level. This ad-hoc definition takes account 
of the fact that the lowest coefficients are the most 
accurately determined and allows errors to be es- 
timated directly. Fig. 3 shows the quark mass de- 
pendence of the onset \x. Note that, from a study 
of the Lee- Yang zeros in the complex m-plane, we 
expect the system to undergo a 'deconfinemcnt 
transition' in the region m < 0.01 on a 6 4 lattice 
at P = 5.1. The scaling with y/m is clearly pre- 
ferred, consistent with the onset of the transition 
being controlled, in part, by a Goldstonc boson. 

2. Conclusions 

Lattice QCD at finite baryon density contin- 
ues to exhibit a phase transition inconsistent with 
He ~ ^ despite the inclusion of the complex 
fermion determinant. It may well be that a suc- 
cessful formulation of QCD at finite density must 
override the effects of a Goldstone boson. 
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Figure 2. The fermion number density 
on 6 4 (triangles) and 8 4 (stars) lattices at 
= 5.1, ma = 0.01. 
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